Stable spatiotemporal solitons in Bessel optical lattices 
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Abstract 

We investigate the existence and stability of three-dimensional (3D) solitons supported by cylin- 
drical Bessel lattices (BLs) in self-focusing media. If the lattice strength exceeds a threshold value, 
we show numerically, and using the variational approximation, that the solitons are stable within 
one or two intervals of values of their norm. In the latter case, the Hamiltonian-vs.-norm diagram 
has a "swallowtail" shape, with three cuspidal points. The model applies to Bose-Einstein con- 
densates (BECs) and to optical media with saturable nonlinearity, suggesting new ways of making 
stable 3D BEC solitons and "light bullets" of an arbitrary size. 

PACS numbers: 42.65.Tg, 42.65.Jx 
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Nonlinear wave dynamics in lattice potentials has well-known realizations in a wide range 
of physical settings, such as localized modes in solid-state lattices jl] and waves in Bose- 
Einstein condensates (BECs) trapped in optical lattices (OLs) Q. Important realizations 
were predicted in optics, in terms of the light transmission in arrays of coupled nonlinear 
waveguides |3[ and optically-induced structures in photorefractive (PR) materials (see sur- 
veys of recent results in Refs. Q and Q), in layered slab media_|6|], etc. Stable solitary pulses 
trapped in lattices in dissipative media were also predicted (7[. Experiments demonstrate 
that lattices in PR crystals indeed give rise to localized quasi-discrete dynamical patterns, 
such as fundamental spatial solitons [sj], vortices [^J, and soliton trains jlflj ]. 

Optical spatiotemporal solitons (alias "light bullets" ) have attracted of lot of attention in 
the last several years, see a review They are 3D wave packets in which the dispersion 
and diffraction are simultaneously balanced by nonlinearity. The search for suitable media 
for the creation of stable "bullets" is a challenging problem, as multidimensional solitons in 
Kerr-type focusing media are unstable against collapse (blowup) Q. The same problem 
impedes creation of multidimensional solitons in self-attractive BECs. Various schemes to 
stabilize solitons in cubic media were proposed. In terms of BECs, a promising avenue is 
the use of periodic OL potentials. 2D and 3D solitons can be stabilized by OLs having the 
same dimension [13j, or by low- dimensional OLs , i.e ., 2D and 3D solitons may be stable in 
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a ID lattice can support stable 3D 



the presence of a ID and 2D lattice, respectively 
solitons in combination with the time-periodic alteration of the nonlinearity sign, provided 
by the Feshbach resonance [lfi| . An essential advantage offered by the low-dimensional OLs 
is that the solitons can freely move in the unconfmed direction(s), which opens a way to 
study their collisions and bound states lj, Q • Essentially new possibilities are opened in 
optics and BECs by radial OLs with axial symmetry, that can be induced by diffraction-free 
cylindrical Bessel beams ^|. Several 2D soliton families, including fundamental solitons 
trapped at the center or in different rings of the Bessel lattice (BL) and dipole-mode soli- 
tons, were predicted in a cubic nonlinear medium, and were shown to be stable in certain 

n 

parameter domains [18j . Very recently, stable ring-shaped solitons in cylindrical lattices 
were also predicted in media with the saturable nonlinearity of the PR type ^| . 

The difference between the previously studied 2D periodic OLs and the BL is more than 
just another symmetry. First, the BL makes it possible to create a soliton at a prescribed 
location (if the experimental field is large enough, several solitons can be created and ma- 
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nipulated simultaneously). Moreover, each soliton can be moved across the field by a slowly 
sliding laser beam. Soliton networks and wires can be also created with arrays of Bessel 
beams, a goal hard to achieve with harmonic lattices. Quite important is also the difference 
in the spatial scales. Indeed, the period of OLs is dictated by the laser-beam's wavelength, 
Ao (although the period can be made larger by launching the interfering waves under an an- 
gle to each other p). On the contrary, the radial scale of the BL is completely independent 
of Ao, being determined by the phase mask through which the beam is passed. The use of 
the BL with an essentially large scale will make it possible to create big paraxial solitons in 
optical media, and big solitons in BEC. The scale difference also highlights the distinction 
of the BL from the single- mode optical fibers, or cigar-shaped traps in BEC (20(, where the 
radius is also limited to microns. 

The subject of this work are 3D solitons trapped in the BL, in the case of the cubic 
self-attraction. This configuration can be directly implemented in BECs with a negative 
scattering length, illuminated by an optical Bessel beam. The corresponding normalized 3D 
Gross-Pitaevskii equation (GPE) for the wave function q is Q 

dq 1 f&q d^q ffiq\ _ , 

where £ is time, r and (77, £) are, respectively, the coordinates along the beam and in the 
transverse plane, r = [rf + C 2 ) 1 ^ 2 , an d the potential function II(r) is proportional to the 
local intensity of the red-detuned optical wave which traps atoms at its anti-nodes (2J [being 
interested in self-supporting solitons, we do not include a magnetic (parabolic) trapping 
potential]. In this work we mainly concentrate on the case of II(r) = vJq (V^Pr), with the 
Bessel function Jo generated by the diffraction- free cylindrical beam |Tj|. Constants p and (5 
determine the strength and radial scale of the BL. Note that the effective potential created 
by the Bessel beam proper is proportional to — J5 2 , while the case of II ~ J actually implies 
interference between the cylindrical beam and a stronger plane wave with an amplitude 
Aq at the same frequency. Indeed, the latter configuration implies the full field amplitude 
A(r) = Aq + A\e tx Jo (^s/2/3r) (with constant A\ and a phase shift x)i which generates 
the above potential through the expansion \A(r)\ 2 w A 2 + 2A A 1 (cosx)Jo (V^r). We 
demonstrate below that the 3D solitons supported by the potentials — J5 2 and — Jo may differ 
because of different degrees of localization of these potentials. The latter potential (which 
was not proposed before) may be more convenient, as it may be controlled by the phase 



shift x) i n particular, it is possible to flip the sign of the Jo potential, which is impossible 
for its Jq counterpart. 

In nonlinear optics, £, (^,C)> an d T are > respectively, normalized propagation and trans- 
verse coordinates and local time ll| , assuming anomalous chromatic dispersion. Then, Eq. 
(fT| governs the nonlinear transmission of a probe signal in a strong BL induced by a pump 
beam launched in an orthogonal polarization, or at a different wavelength. In this case, 
Il(r) measures local modulation of the refractive index for the probe signal induced by the 
pump beam through the cross-phase modulation, and the choice of II ~ J again implies 
the interference between the Bessel beam and a stronger plane wave. The pump beam is 
assumed to propagate in a nearly linear regime, that actually implies strong saturation of 
the pump beam in the medium with saturable nonlinearity. 

Equation conserves the norm, U — J J f \q(ri,(,T~)\ 2 dr/d^dr, Hamiltonian, H, and 
angular momentum. We search for stationary solitons, q(r], (, r, £) = w(r, r) exp(i6£), where 
b is a real propagation constant, and w, which must vanish at \t\,t — > oo, obeys the equation 

w TT + w rr + T~ x w r + 2[pII(r) - b]w + 2w 3 = 0. (2) 

In numerical calculations, we fix the radial scale, (3 = 10, and vary b and p (cf. Ref. [rs|). 

Using an obvious Lagrangian representation of Eq. (J2J), we start with the variational ap- 
proximation (VA). To this end, we adopt the Gaussian ansatz for the 3D soliton, w(r,r) = 
Aexp [— (r 2 /T 2 + r 2 /R 2 ) /2], with free parameters T, R and A. The known procedure 



2l| leads to relations T 2 = R 2 / (l - ppR 4 e~^ 2 / 2 ), b = p(l- 3/3R 2 /A) e"^ 2 / 2 + (4R 2 ) \ 
and U 2 = 8vr 3 i? 2 (l - pf3R 4 e~P R2 / 2 ), where U = ti 3 / 2 A 2 TR 2 is the soliton's norm. These 
equations predict the main characteristic of the soliton family, U(b). Then, the Vakhitov- 
Kolokolov (VK) criterion, dU/db > 0, provides for stability of the solitons against pertur- 



12 1 . A first VK-stable interval appears around R = R { 



(i) 



bations with real eigenvalues 
2.75/v^, when the BL strength attains a threshold value, p^. w 0.12/?. A second stability 

1/2 

interval appears around i% = 2 ^3 — v^3") /(3 ~ 1.59/v 7 /?, when p passes a second 
threshold, p\^ r ~ 0.32/5. As shown below, numerical results confirm the existence of two 

(2) (1 2) 

stability intervals for p > p t ^ r . It is relevant to compare the above values Rq ' and the size 
of the central core of the BL potential in Eq. (J2J), which is determined by the first zero of 
the function J (^/2f3r^, r w 1. 70/>//7. The comparison suggests that the above simple 
ansatz, that neglects BL-induced fringes in the shape of the soliton, is appropriate, as the 
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predicted solitons do not cover many rings of the radial lattice. 

Numerically, soliton families were found from Eqs. (J2J) by means of the relaxation method. 
Simulations of the full GPE (|TJ were carried out using the Crank-Nicholson scheme. To 
test the stability, we took perturbed solutions as q = [w(r, r) + u(r, r) exp(5£ + inifi) + 
v*(r, t) exp(<5*£ — in4>)\ exp(z&£), where <fi is the azimuthal angle in the plane (77, £), and 
n = 0,1,2,... . The linearization of Eq. (JIJ) leads to straightforward equations for the 
perturbation mode (u, v) and its growth rate 5. The numerical results are summarized in 
Fig. 1. We found that the 3D solitons exist in the Jo potential if the propagation constant b 
exceeds a certain cutoff value b co . Figure 1(a) shows that b co increases with the BL strength 
P [bcoip) vanishes at p = as unstable free-space 3D solitons exist for all b > 0]. The 
soliton's norm U is a non-monotonous function of b, see Fig. 1(b). The prediction of the 
VK criterion, that parts of the soliton family with dU/db > may be stable, was verified 
by numerical computation of the growth rates, Re(5). It was found that the solitons are 
completely stable precisely in regions singled out by the condition dU/db > 0, see Fig. 1(c). 
For the Jo potential and relatively small values of the BL strength p, the solitons are stable 
in a narrow interval of b to the right of b co [see a small solid segment of the curve in Fig. 
1(b) corresponding top — 5]. Losing their stability for b exceeding the value at which dU/db 
vanishes, the 3D solitons drastically differ from their BL-supported 2D counterparts, that 
may be stable in their entire existence domain if the lattice strength p is large enough 3| • 
Remarkably, for moderate values of p, the 3D solitons are stable in two disjoint intervals of 
b [one abuts on the cutoff point, b = b co , and the other is found at relatively large values of 
b, see the curve corresponding to p — 5.5 in Fig. 1(b)]. This feature was predicted above by 
the VA. However, there is only one soliton-stability domain, adjacent to the cutoff, for the 
J5 2 potential, see Fig. 1(d). The difference can be traced to different localization degrees of 
the Jq and J potentials. 

The Hamiltonian-vs.-norm diagrams for the soliton families, which is a useful tool for 



the analysis of the soliton stability 



22j, are plotted in Fig. 2. With the Jq potential, they 



exhibit one or three cuspidal points. In the latter case [Fig. 2(a)], the diagram features a 
"swallowtail" pattern, which accounts for the existence of the two above-mentioned distinct 
stability regions for the 3D solitons. Although this pattern is one of generic possibilities 
known in the catastrophe theory, it rarely occurs in physical models [22]. In Fig. 3 we 
display three typical examples of stable 3D solitons supported by the Jq potential. A low- 



5 



amplitude broad soliton covering several lattice rings, with the propagation constant b close 
to the cutoff b co , is shown in Fig. 3(a), and high-amplitude narrow solitons, mostly trapped 
within the BL core, are presented in Figs. 3(b,c). 

An important issue for the 3D solitons in the focusing medium is the occurrence of 
collapse. We expect that the solitons which are unstable against small perturbations either 
undergo the collapse, or decay into radiation waves, depending on the perturbation. To 
test these expectations, we simulated the evolution of perturbed solitons, starting with 
Q{€ = 0) = w(r], £, t)(1 + ef), where e is a small perturbation amplitude, and / is either 
a random variable uniformly distributed in the interval [—1,1], or constant, f — 1. We 
have found that the solitons which were predicted to be stable in terms of the perturbation 
eigenvalues are indeed stable in the direct simulations. As an example, in Fig. 4 we show the 
evolution of a soliton with a random perturbation whose amplitude is e = 0.05. Despite a 
considerable transient change of its shape, the soliton relaxes to a self-cleaned form which is 
quite close to the unperturbed one. Under the same random perturbation, linearly unstable 
solitons were found to decay into radiation. The uniform initial perturbation, / = 1, amounts 
to an instantaneous change of the soliton's amplitude. This perturbation excites persistent 
oscillations of linearly stable solitons, while the effect on linearly unstable ones is completely 
different from the response to random perturbations: they do not decay, but, instead, blow 
up. 

In conclusion, we have constructed families of 3D solitons supported by zero-th order 
Bessel optical lattices. We have found that, in one or two intervals, the solitons are stable, 
provided that the lattice is strong enough. We also expect that 3D solitons can be stabilized 
by higher-order Bessel lattices, too. This is the first example of a 3D soliton stabilized by 
an axially symmetric lattice; unlike stable 3D solitons in 2D periodic lattices, the soliton's 
size is not determined by the wavelength of the laser beam, but may be chosen at will. 
Depending on the lattice strength, the Hamiltonian-vs.-norm diagram for solitons in the Jo 
potential displays a "swallowtail" bifurcation, a rare phenomenon in physics. 

The results suggest new approaches to the still unresolved 111] problem of the creation of 
3D solitons in the experiment. One possibility is to make a BEC soliton in the cylindrical 
optical lattice, which may be controlled by the superposition of the Bessel beam and uniform 
background field. In optics, creation of the 3D lattice solitons require the bimodal setting in 
a medium with saturable nonlinearity. In all these settings, stable 3D solitons are predicted 
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to feature universal properties. 
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FIG. 1: (a) The propagation constant cutoff vs. the lattice strength, shown for the cylindrical 
potentials of both types, Jo and J5 2 . (b,d) The soliton's norm vs. the propagation constant for the 
potentials Jo and J5 2 , respectively, (c) The perturbation growth rate vs. the propagation constant 
for the Jo potential. Here and below, solid and dashed curves correspond, respectively, to stable 
and unstable branches. 

FIG. 2: The Hamiltonian-vs.-norm diagrams for p = 5.5 (a) and 8 (b). Here and below, all the 
figures are shown for the Jo potential. 

FIG. 4: The evolution of a stable soliton under a random perturbation, for p = 5, b = 0.11, 
U = 6.477: (a) £ = 0, (b) £ = 200, (c) f = 500. 



FIG. 3: Isosurface plots for stable solitons: (a) p = 5, b = 0.11, U = 6.477; (b) p = 5.5, b = 0.85, 
U = 3.677; (c) p = 6, b = 2, U = 3.591. 
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